Depletion forces near a soft surface 
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We investigate excluded-volume effects in a bidisperse colloidal suspension near a flexible interface. 
Inspired by a recent experiment by Dinsmore et al. (Phys. Rev, Lett. 80, 409 (1998)), we study 
the adsorption of a mesoscopic bead on the surface and show that depletion forces could in principle 
lead to particle encapsulation. We then consider the effect of surface fluctuations on the depletion 
potential itself and construct the density profile of a polymer solution near a soft interface. Sur- 
prisingly we find that the chains accumulate at the wall, whereas the density displays a deficit of 
particles at distances larger than the surface roughness. This non-monotonic behavior demonstrates 
that surface fluctuations can have major repercusions on the properties of a colloidal solution. On 
average, the additional contribution to the Gibbs adsorbance is negative. The amplitude of the 
depletion potential between a mesoscopic bead and the surface increases accordingly. 



I. INTRODUCTION 

In a recent experiment, Dinsmore and collaborators ^ 
measured the probability distribution of a single colloidal 
bead enclosed in a rigid, pear-shaped vesicle. They find 
that upon addition of much smaller spheres, density gra- 
dients give rise to an effective force field that pushes the 
particle towards the wall. These depletion forces have 
been discussed for almost fifty years I^IE^I^, especially 
since they can be used to organize self-assembled struc- 
tures The remarkable feature observed in T] is that 
entropic interactions drastically affect the distribution in 
a way that depends on the local shape of the membrane. 
Indeed, it has been shown ^ that inhomogeneities in 
the curvature distribution induce additional forces that 
bring the particle to inwardly curved regions. The au- 
thors then suggested that if curvature induces forces on 
particles, the latter may in turn lead to shape transitions 
provided that the membrane is flexible enough Q]. 

If this last statement is valid, it would be conceivable 
to take advantage of depletion forces to mimic biological 
processes such as protein integration, membrane fusion, 
or encapsulation. For instance, polyethylene glycol is a 
(polymeric) depleting agent known for its low toxicity: 
By changing its mass or concentration, one can selectively 
tune the range and the depth of the depletion potential 
without affecting any physiological parameter, i.e., pH, 
temperature and ionic strengths. This remarkable ver- 
satility further justifies the growing interest in depletion 
interactions and makes them particularly suitable for ap- 
plications such as gene transfection 8]. 

In this article, we present a systematic study of en- 
tropic forces near a fluctuating interface. The additional 
degree of freedom associated with surface softness in- 
duces numerous antagonistic features that have to be 
carefully accounted for. Besides the eventuality of shape 
transformation, a fluid membrane experiences thermal 



undulations whose amplitude is inversely proportional to 
its bending rigidity j^. The steric forces between adja- 
cent bilayers, originating from shape excitations, are in- 
voked to explain the phase behavior of lamellar stacks of 
membranes ^|. Similarly, a colloidal particle has first 
to overcome a repulsive barrier prior to adhering to a 
fluctuating membrane. This point is usually ignored in 
discussions concerning colloidal adsorption, even though 
it may deeply affect the kinetics of the binding process. A 
more refined consequence is the effect of membrane fluc- 
tuations on the concentration profile of small particles, 
which in turn affects the depletion potential. It is the 
aim of this work to theoretically address these different 
issues. The paper is organized as follows: In Section ITTI 
we give a brief outlook of some general results regard- 
ing effective, entropic interactions. The balance between 
depletion and bending forces is studied in Section IIIII 
where we derive a tentative phase diagram of the sys- 
tem. We present in Section IIVI our findings concerning 
the corrections to the depletion potential. We conclude 
this report with a discussion and summary of our results 
m Section El 



II. EFFECTIVE INTERACTIONS IN A 
COLLOIDAL SOLUTION 



We consider a binary mixture of small and large spher- 
ical beads with radius Rg and R, respectively. In a single- 
component description, the large particles can then be 
seen as movi ng i n an effective force field created by the 
other species From a theoretical point of view, this 
corresponds to integrate out the partition function over 
the degrees of freedom of the small particles. The par- 
tial trace can be performed for a dilute solution within 
the formalism of density functional theory low 
density, the effective potential W{v) acting on a single 
colloid immersed in a "sea" of small particles is [T^ 
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with (3 = {kBT)~^, p(r) the density of small compo- 
nents and pi, — p(oo) the bulk value. The Mayer function 
equals /(r) = exp[— /3u(r)] — 1, with i;(r) the direct in- 
teraction potential between small and large components. 
Note that in Eq. , the thermodynamic quantities have 
to be evaluated before insertion of the large particle. This 
assertion happens to be essential since the symmetries of 
the relevant functions are solely determined by the ge- 
ometry of the system. For instance in the presence of a 
planar wall at z = 0, the density profiles depend on the 
z-coordinate only, whereas translational symmetry in the 
{x, y)-plane is broken after particle insertion. 

Eq. (Q) provides a low-density expansion for the deple- 
tion potential even if the interactions between the species 
are soft and possibly contains attractive contributions. It 
also encompasses the hard-sphere limit considered in the 
original theory of Asakura and Oosawa , for which the 
Mayer function reduces to /irs(r) = —H{R + Rg — |r|). 
H(x) is the usual Heaviside step function that equals 
1 for X > and otherwise. In the presence of a 
planar wall, the volume fraction of depleting particles 
(/7(z)=47ri?^p(z)/3 is also expressed in terms of the step 
function, namely ip{z) = ipi,H(z — R^). The convolution 
Eq. Q is then a measure of the overlap volume — see 
Fig. 1. For a separation I between the surface of the large 
sphere and the wall, the depletion potential reads P.[l^ 



/3M^(o)(?) = -Sip, 




for I < 2Rs, and = otherwise. The value 

at contact is /3Vl^(")(0) = ~3ipta-\ where we introduce 
the aspect ratio a — Rs/R- For typical volume fraction 
(fib — 0.2 and aspect ratio a — 0.1, the depth of the po- 
tential well is VF(°)(0) = -eksT. Thus, these entropic 
effects cannot be ignored as they may eventually over- 
come the direct interactions. Actually, they appear to be 
essential in describing phase separation and aggregation 
phenomena in polydisperse mixtures . 

While much of the theoretical and numerical effort as- 
sumes that the small particles are structureless spheres, 
experimental studies mainly use polymers as depleting 
agents. Actually, the formal mapping of a polymer solu- 
tion onto a fluid of "soft" particles is a challenging prob- 
lem because of intricate many-body interactions [la . Il7j . 
For most of practical purposes however, the exact na- 
ture of the small components can be ignored: at low 
concentrations, polymer coils can be regarded as spheri- 
cal particles of radius Rg, their radius of gyration. Note 
that this assertion is still valid above the overlap con- 
centration, provided that one substitutes Rg with the 
correlation length of the solution . 

In the recent years, a variety of experimental tech- 
niques have been developed to probe soft materials at the 
sub-micrometer scale: optical tweezers, surface force ap- 
paratus (SEA), atomic force microscopy (AFM), or total 
internal reflection microscopy (TIRM). These progress 
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FIG. 1: Illustration of the depletion attraction phenomenon. 
The center-of-mass of the small particles is excluded from a 
shell of thickness Rs around the large spheres and along the 
surface (delimited by the doted lines). When these depletion 
volumes overlap (grey regions), the osmotic pressure of the gas 
of small particles becomes anisotropic and the large particle 
is pushed toward the surface. The effective volume available 
for the small particles is further increased if the surface can 
deform and envelop the large particle. 

have enabled direct comparison with theoretical pre- 
dictions, and numerous experimental works have con- 
firmed the validity of Eq. jSJ for volume fractions (pb < 
0.25 0,0, 1131 • At higher concentrations, a repulsive 
barrier forms at distances I ~ 2i?s, followed by secondary 
minima reminiscent of the la yer ed structure of the small- 
component hquid mmiUlilllil. Eventually, kinetics 
effects will prevent the colloidal particle to reach the pri- 
mary minimum In the present work, we assume 
that the packing fraction is low enough to neglect the 
correlations between small particles. The kinetics is also 
ignored throughout the treatment p^. 

III. DEPLETION-INDUCED WRAPPING 

The adsorption of colloidal particles on a soft surface 
is an essential step in many biological processes. The un- 
derlying physics of this mechanism has been extensively 
studied for simple model systems [2^ ll^, 113, 1^ ■ Here, 
we investigate the special case of adhesion driven by de- 
pletion forces. We specify the interface profile through 
the displacement field u(r||), where r|| = ix,y) spans the 
horizontal reference plane at z = 0, and define the Fourier 
representation Uq = / d^ry exp(— iqrj|)w(r||). The sta- 
tistical behavior of the surface can be studied within a 
continuum model based on the following Hamiltonian 

Wx=i/ {r + aq' + .q') , (3) 

with K the bending rigidity and a the surface tension. 
The elastic coefficient r is related to the curvature of an 
eventual external potential (gravitational field, effective 
"mean- field" intermembrane interactions). This model 
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allows us to describe a wide class of surfaces, ranging 
from liquid-liquid or liquid-gas interfaces, with high sur- 
face tension and vanishing bending modulus, to pure fluid 
membranes with finite bending modulus and zero surface 
tension. The intermediate combination where both a and 
K are finite would for instance correspond to a surfactant 
monolayer at the oil/ water interface !32]. 

A flexible surface as described by the Hamiltonian T-Cj 
is not a simple planar surface, but i) can be deformed 
under an external constraint, and ii) presents an intrinsic 
roughness due to thermal fluctuations. In this section, 
we derive the phase diagram for the different wrapping 
states of the colloid. For the sake of simplicity, we restrict 
ourselves to the bending regime r = = 0. The critical 
behavior of the system results from a balance between 
the driving depletion force and the membrane bending 
energy. The former separates into an entropic and an 
enthalpic contribution, that we study separately. 



A. Depletion vs. Helfrich repulsion 

A colloidal particle feels in the vicinity of the mem- 
brane a steric repulsion resulting from the thermal undu- 
lations All analytical expression for this repulsion 
is known only for a membrane close to an infinite, flat 
wall: the free energy at given mean separation / scales 
as VHelfr^chil) = AikBT)y{nP), with A = 7rVl28 [13. 
One way to make contact with the curved geometry is 
to use the Derjaguin approximation |3^ . which relates 
the force between a sphere and a plane to the poten- 
tial between two planes. In the limit I <^ R, the force 
exerted by the fluctuating surface on the spherical par- 
ticle is F{1) = 2'KRVHeifrich{l)- Here, I is the dis- 
tance between the surface of the particle and the aver- 
age position of the membrane. It follows that the effec- 
tive potential between a sphere and a membrane equals 
Vrep{l) = /;°° F{z)dz. For / < 2Rs < i?, the total poten- 



tial V{1) =Vrep{l) + W^^\l) 
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with k = K./{kBT), whereas only the repulsive part sur- 
vives at separations I > 2Rs. Minimizing V gives the po- 
sition of the minimum « (7r^/(192 More- 
over, the requirement that l^, corresponds to a global min- 
imum defines a critical bending modulus 



(5) 



■ 192 (^h 

Below this critical value, there is no bound state and 
the particle can never reach the surface. In the adhesive 
regime k > kc,i, the particle localizes at the minimum l^ 
of the potential that equals le ~ 0.5i?s for k = 10 and 
ipb = 0.2. 

One could argue that the superposition of the Helfrich 
repulsion and the direct potential in an attempt to pro- 
duce an effective potential is in general not correct [35|- 



Rather, a virial expansion of the direct part is needed to 
describe the thermodynamics of the system [s^. Here 
however, both the attractive and repulsive parts of the 
potential are entropic in nature, thus possibly justifying 
such a superposition. 



B. Depletion vs. bending energy 

We now assume that once the repulsive barrier is over- 
come, the colloidal particle is adsorbed on a purely elastic 
sheet. As shown on Fig. 1, the depletion contribution to 
the free energy is further decreased if the membrane en- 
velops the particle. Note that we assume here that the 
solvent on both sides of the surface is incompressible: the 
only volume contribution comes from the variation of the 
overlap volume. The natural parameter that describes 
the wrapping transition is the degree of penetration lo 
defined as the ratio between the covered area and the to- 
tal area of the colloid. The equilibrium value of uj is set 
by the balance between osmotic and bending forces. Ne- 
glecting the fluctuations, the free energy contributions 
divides into three parts. The first term is the bending 
energy of the covered area that simply equals Sttkco. In 
the limit of a small aspect ratio, the overlap volume is 
^Vov = SirR'^LLia + O(a^), so that the depletion- induced 
adsorption can be described in terms of a constant ad- 
hesion energy per unit area F = 2kBTptRs This 
approximation is certainly not valid in the early stages 
of the wrapping process where lo a, but remember that 
we are interested by the limit w ~ 1. The last contribu- 
tion arises from the free part of the membrane and is 
more difficult to evaluate, since it generally requires one 
to solve the full non-linear shape equations ^3 US ■ 
Fortunately the problem greatly simplifies for a mem- 
brane with no surface tension, providing that there is no 
pressure difference between both sides of the surface. In 
this particular case, no forces are applied to the free part 
of the membrane. The deformation profile becomes a 
catenoid, a surface with zero curvature and consequently 
a vanishing bending energy. The free energy of the sys- 
tem is then a linear function of w 

(3AF{w) = Snuj (^k - ^ipba-^^ . (6) 

Eq. 10 then defines a second critical bending rigidity 

kc,2 = -^ipta"'^ . (7) 
An 

The colloidal bead is in the unwrapped state above Kc.2, 
whereas the membrane completely encapsulates the par- 
ticle below this critical value [i^. It is worthwhile to 
note that the minimum of the free energy is not an ex- 
tremum, but corresponds to either boundaries of the do- 
main Lu = or 1. Therefore, the equilibrium profile 
does not have to satisfy the usual stationarity require- 
ment c* = l/i?+ ^2r/^ 
the detachment point |37l | 



with c* the mean curvature at 
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FIG. 2: Schematic phase diagram for the wrapping transition. 
(U): unbound state. (B): bound state. (W): wrapped state. 



The conditions Eq. ||SJ and |7J) enable us to draw the 
phase diagram of the system presented in Fig.|5| For low 
bending rigidity or volume fraction, the colloid is always 
in the unbound state because of the Helfrich repulsion. 
Also, the membrane cannot wrap the particle if the bend- 
ing modulus is too high. Conversely, the membrane en- 
velops the colloid when the bending modulus lies within 
the range Kc.i < H < Kc,2- For typical values a = 0.1 
and iph = 0.2 the boundaries correspond to Kc,i ~ 0.77 
and Kc,2 ~ 4.77. Depending on the volume fraction and 
the aspect ratio of the small colloids, this interval can 
actually vanish: the condition of existence of a wrapped 
state is given by the inequality 



an 



(8) 



Arguably, the heuristic derivation presented here is not 
accurate enough to predict the nature of the transitions 
and the exact values of the numerical constants. The in- 
fluence of finite-size effects on the unbinding transition 
would in particular deserve a further analysis in the fu- 
ture. However, this phenomenological description pro- 
vides some important insights for the understanding of 
the wrapping transition. 



IV. FLUCTUATIONS OF THE INTERFACE 
AND CORRECTIONS TO THE DEPLETION 
INTERACTION 

We derived in the preceding section the phase diagram 
for the wrapping of a colloidal particle. In particular, 
we found that encapsulation can indeed be driven by de- 
pletion forces. However, we have not yet addressed an- 
other relevant question arising from the surface softness, 
namely the effect of thermal undulations on the depletion 



potential itself. As explained in Section H, the deriva- 
tion of this potential only involves the density profile of 
the small components. Here we adopt a procedure in a 
sense opposite to the original route pursued by Asakura 
and Oosawa, and construct the equilibrium concentration 
profile for a polymer solution, instead of spherical beads, 
near a fluctuating surface. We then apply this result to 
evaluate the effective depletion potential. 



A. Statistical mechanics 

We consider a generic interface described by the Hel- 
frich Hamiltonian Eq. Q. The thermodynamic prop- 
erties are specified by two lengths ^|| and [Sg: ^|| 
characterizes the exponential decay of the height-height 
correlation function 



G(r||) = (^r||HO))^^ e 



(9) 



whereas — ( u(r||)^)^''" is a measure of the roughness 



of the surface. For instance, a fluid membrane has a 
vanishing surface tension a — and the in and out-of- 
plane correlations lengths are ^|| = (4K/r)^/'* and = 

^I|/(4a/k), respectively. 

If a colloidal solution is in contact with this surface, 
we argue that the rigorous treatment of the fluctuations 
should be as follows: assume that the density profile p[u] 
is known for a given conformation M(r||) of the surface. 
The thermal fluctuations can then be traced out accord- 
ing to 



p{z) = (/7(r||,z))^ 



Vue 



-m-i 



p[u\ 



(10) 



so that p{z) is finally expressed as a functional of the 
Fourier transform of the surface correlation function 



G(q) - (Kl'h 



(11) 



The remaining problem is to evaluate the functional 
p[u\. To this end, we will assume that the depleting 
particles are linear polymer chains instead of spherical 
beads. This idea might at first seem questionable as a 
macromolecule has many more degrees of freedom than 
a simple hard sphere. Nevertheless, the partition func- 
tion 2^jv(r|| , z) of an ideal chain of N monomers with one 
extremity fixed at (r|| , z) can be obtained using an anal- 
ogy to Brownian motion |4ll |. The complexity of surface 
impenetrability then reduces to a boundary condition 
problem that can be solved analytically. The polymer 
concentration is then expressed in the low-density limit 
as El 



N-n{-r\\,z) 



(12) 
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Far from the wall, the bulk value equals pb = Af /V with 
M the total number of chains in solution. The average 
value Eq. H10|l of the density profile follows immediately. 

The evaluation of Zj^ for a an y su rface profile u(r||) is 
in general a fairly complex task psLli^ . Here, we follow 
a perturbative approach based on the work of polymers 
grafted to a fluid membrane published elsewhere pslli^ . 
The key ingredient of the derivation is to write the par- 
tition fmiction as a series -Eat = X^i^^^j where z''^!^ is 
of order u'. The expansion of the partition function of 
a Gaussian chain in the vicinity of a non-planar surface 
is detailed in the Appendix. All the results displayed in 
the following are exact up to linear order in the natural 
expansion parameter k^^ . 

At this point, it should be mentioned that our pertur- 
bative scheme is only valid in the limit of "almost flat 
surface", i.e., for much smaller than the radius of gy- 
ration Rg of the coils. This condition can be written in 
terms of the in-plane correlation length as 



(13) 



It implies that our calculations are valid even for Rg < ^|| , 
provided that the bending rigidity is large enough. For 
instance, for a bending modulus k = 10, the consistency 
of the approach is still fulfilled for ^|| as high as ^|| = 
10i?„. 



B. Gibbs adsorption 

We begin the discussion with the first moment of the 
distribution, which is less computationally demanding. 
In the vicinity of the interface, the proportion of polymers 
that belongs to the adsorption or depletion layer may be 
quantified by the adsorbance F defined as |43 



dz ipiz) - Ply) 



(14) 



The adsorbance usually scales as F ^ ^pf,, with ^ the 
typical length-scale over which the concentration relaxes 
towards the bulk value. Note that the length ^ is an 
algebraic quantity: a positive value reflects particle ac- 
cumulation, whereas a negative value expresses a deficit 
of particles in a layer of thickness ^. For instance, the 
concentration profile p^"-* of ideal chains depleted from a 
flat, infinitely repulsive surface is [iMIi^ 

P^'\0 =Pb[e [l + erf(C/2) - 2erf(C)] 2erf(C/2) 



erf(C) - -r 
a/tt 



(15) 



with C = "^/Rg the (dimensionless) distance from the 
wall, and Rg — \/ Na^ /Q the radius of gyration of the 



polymer coil. The corresponding adsorbance is [5C 

2 



r(") 



-PbR. 



a ■ 



(16) 



In the rigid wall limit {k — > oo), the particles are ex- 
cluded from the vicinity of the surface over a distance 
doo ~ Rg- If the prescription R^ = 2Rg/^ is used, 
the excess surface density for the hard sphere and the 
polymer solutions are the identical, which illustrates the 
connection with the original model of Asakura and Oo- 
sawa. 

For a finite rigidity (k < oo), one would at first expect 
the thickness of the depletion layer to increase according 
to da ~ Rg + £,^_■ This can be checked by performing the 
trace over the surface conformations Eq. I|10(l . We find 
the correction AF = F — F*^") to the adsorbance 



AF 



1 Pb 

^/t^ Rg 



{2n) 



-I- ^/TrqRgCTf{qRg) - 1 



(17) 



where the correlation function G(q) is defined by 
Eq. One should note that AF is always negative, 

reflecting an expansion of the depletion layer. Expres- 
sion (|17|l can be evaluated for a membrane without sur- 
face tension: we get for cr = 



AF 



-K^^PbRg 



for < Rg, 
for > Rg. 



(18) 



We see that in both limits, the depletion layer is increased 
by min(^|| , i?g). This seems to confirm the "naive" 
expectation that arises from the scaling argument. Nev- 
ertheless, this picture happens not to be completely ac- 
curate, since the adsorbance is only the first moment of 
the distribution and therefore does not contains all the 
information about the density proflle. 



C. Correction to the density profile 



More difficult to evaluate analytically is the correction 
Ap{z) = p{z) — p^^\z) to the concentration proflle. We 
flrst focus on the proximal region z <C Rg [5l| . The value 
at the origin is easily obtained and equals 



Ap(0) = ^ (f 

16k \Rg 



(19) 



Quite unexpectedly, we flnd that the polymers accumu- 
late near the surface: Ap(0) > 0. Moreover, if we extend 
the calculations up to linear order in z we obtain for a 
generic surface 



+ AV^qRgerfiqRg) - l] 



(20) 
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FIG. 3: Correction to the density profile for an interface with 
vanishing surface tension ((7 = 0). The different curves corre- 
spond to C||/^s = 0.1, 0.5, 1 and 2, respectively. In the proxi- 
mal region, membrane fluctuations favor polymer adsorption 
on the surface; at larger distances, the chains are pushed back 
into the bulk by the Helfrich repulsion. 



Since the integrand is positive for any value of the param- 
eters, the correction to the depletion profile starts with a 
negative slope. For a membrane without surface tension, 
we get the following limiting cases 



Ap(z) 




2%/2z 



for Rg > 
for Rg < f II 



(21) 



This shows that the simple picture emerging from the 
scaling argument is in fact erroneous. Formally, the be- 
havior described by Eqs. (|19|l - (|21|l is reminiscent to the 
adsorption of ideal polymer chains on a flat wall ■ It 
should be emphasized that the effects discussed here are 
more than a simple correction to p^*'^ Indeed, expanding 
Eq. (|15|l for z <^ Rg gives p*^°)(z) oc z^, which vanishes 
at the surface. Therefore, Eq. (|20|l specifies not only the 
short-distance behavior of the correction Ap but actually 
describes the total concentration profile p — p'"-* + Ap. It 
is even possible to achieve a density at z = comparable 
to the bulk density, but for values of the parameters that 
reach the limit of validity of our approach — see Eq. (|13|l . 

The analytical expression for Ap is uninstructive and 
will not be given here. Fig. 13 presents the correction 
to the density profile for different values of the ratio 
,52iJ. As previously noted, the coils are first ad- 
sorbed in the proximal region, i.e., the chains have the 
tendency to "fill" the thermal undulations. Interestingly, 
the chains are repelled more from the fluctuating surface 
at intermediate distances. This effect ultimately dom- 
inates the overall profile since the adsorbance Eq. H17|) 
is negative. Finally, the interface affects the distribution 
over a finite extension ^ 4i?g. As illustrated in Fig.|21 the 
effect is rather minor for small values of the ratio S^\\/Rg. 
On the other hand, the correction cannot be ignored for 
^11 > Rg, and we plot in Fig. ^ the full depletion profile 
p{z) = p'^^^z) + Ap(z) for a bending rigidity k = 10. 



FIG. 4: Complete density profile p = p'"-* -I- Ap for an interface 
with (T = and k = 10. The different curves correspond to 
(,\\/ Rg = 1, 5 and 10, respectively. 



D. Effective depletion potential 

The non-monotonic behavior revealed in the preceding 
section strongly affects the depletion potential between 
the colloidal particle and the surface, mediated by the 
small polymer coils. At this point of the derivation, we 
have to evaluate the convolution product Eq. (^) numer- 
ically since it involves the complete density profile. For 
the sake of simplicity, we assume that the Mayer function 
characterizing the polymer/colloid interactions is still a 
step function. The correction ISW — W — ly^") is de- 
picted in Fig. 13 in reduced units, where the combination 
PbRg is analogous to the volume fraction cpi, of spherical 

V = 0-2, 



particles. For reasonable values R — lORg, pbR 
K = 10 and = 5Rg, the depth of the well is increased 
by a substantial amount AIF(O) w —lAItBT. In terms 
of the Asakura-Oosawa potential Eq. ||2Jl, this represents 
25% of the original value at the wall. Of course, one has 
certainly to take into account the direct interaction be- 
tween the colloidal particle and the surface at distances 
smaller than , but this question is beyond the scope of 
this report. For ^|| = 5Rg and k = 10, this distance is 
approximatively £,± w OARg. Note also that the curves 
present an infiection point at a separation matching the 
minimum of Ap{z). As increases, this point eventually 
becomes an minimum of AW{z). 

Actually, the amplitude of the effective potential can 
be inferred from very general arguments. Let us start 
with a solution of particles confined between two planes. 
It is easy to show "SS^ that the value w{0) (per unit 
area) of the depletion potential at contact scales as 
(3w{0) ~ r^*'^ with the adsorbance given by Eq. H16|l . 
In the Derjaguin approximation, the corresponding force 
between a spherical particle of radius R and a plane is 
F ~ 'w{Q)R. Since the range of the depletion interaction 
is fixed by the size Rs of the depleting particles, the value 
at contact of the depletion potential between a colloidal 
particle and a plane is PW'^^^O) ^ RRsT'^°\ in perfect 
agreement with Eq. |(2J). Assuming the same scaling for 
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the polymers are closely bound to the surface by an ef- 
fective, attractive potential. This primarily results from 
the tendency of the chains to fill the holes or valleys on 
a rough surface, and can be understood as follows: close 
to a wall, the colloidal solution is affected over a charac- 
teristic length fixed by the size Rs of the particles. This 
naturally gives rise to a surface tension, defined as the 
free energy penalty per unit area to create the depletion 
layer. If the wall is slightly curved with radius R, the 
excess surface energy 7^ can be expanded in powers of 
the ratio q = Rg/R 



(23) 



FIG. 5: Correction to the depletion potential in dimension- 
less units: Afi = {k/ptRl)AW, for ^n/Rg = 1,2 and 5 re- with 7^°) the planar wall contribution that scales as 



spectively. The inflection point is a manifestation of the non- 
monotonic behavior of Ap{z). 



the correction to the potential, we expect AW to be- 
have like PAW{0) ~ RRgAT. It is more instructive to 
consider the ratio AW/W'^^'> 
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(0) 



AM^(O) 
IFW(O) 



[k^^ for ~ Rg, 



(22) 



showing that the effects discussed here can be more than 
a minor correction. 



V. DISCUSSION 

The simple analysis described in this report has pro- 
vided some very interesting and somewhat unexpected 
results for depletion effects at fluctuating interfaces. The 
phenomenological approach initiated in the first part of 
this study reveals that entropic interactions can indeed 
lead to full particle coverage. This route to particle 
encapsulation is particularly attractive since it benefits 
from the wide versatility offered by depletion forces. 

The main part of this work is dedicated to the con- 
centration profile of small particles near a fluctuating 
surface. For the last decade, intensive studies have fo- 
cused on the way a colloidal suspension affects the elastic 
properties of the surface. Particularly well documented 
is the adsorption of polymers from dilute or semi-dilute 
solutions [53, [53, |5g , as well as the depletion of diverse 
objects like r igid rods [sJl, ideal or self-avoiding poly- 
mers 0, [5^. To the best of our knowledge, it is 
the first time that the opposite question, namely, the ef- 
fect of the surface on the bulk properties, is theoretically 
addressed. For this derivation we make three major hy- 
pothesis: i) the particles are ideal polymer chains, ii) we 
treat only the very dilute case where the polymers act in- 
dependently, and iii) we assume that surface undulations 
are only a perturbation to the planar configuration. Even 
within this restrictive set of assumptions, we find the fol- 
lowing remarkable results. In the first regime z < 



ksTpbRs- Note that with this convention, a 
sphere has a positive curvature. The first correction a 
controls the affinity of the solution for the curved surface. 
It can be shown for ideal polymers that a > 0, mean- 
ing that the surface energy is lowered when the chains 
are in the vicinity of a concave region [R < 0). In the 
language of curvature energy, one would infer that the 
surface bends more favorably toward the polymers with 
a spontaneous radius of curvature Rq ~ —Rg- In the re- 
verse problem considered here we can conclude that the 
chains have a higher affinity for the concave regions, re- 
sulting in a finite concentration at z = after averaging 
over the surface degrees of freedom. This effect becomes 
more pronounced as one increases the in-plane correla- 
tion lengths. Indeed, the polymer coils can never fit into 
the fluctuations when <C Rg and the correction is neg- 
ligible. This is not the case in the regime where ^|| > Rg. 
Here, the statistical weight of the undulations with cur- 
vature at their minima smaller than 1 / Rg grows rapidly, 
followed by an increase of Ap(0) as shown in Fig. O 

In the range £,± < z < 4i?g, the chains are on the 
contrary repelled more from the fluctuating surface than 
from a flat wall. This unusual behavior reflects a broad- 
ening of the interface: as the polymer coils move away 
from the average elevation z = 0, they "feel" an effective 
layer of thickness from which they are excluded. This 
repulsion actually corresponds to the force exerted by the 
fluctuating surface on a particle, and suggests that finite- 
size effects may lead to drastic changes to the Helfrich 
repulsion. This question is part of ongoing investigations 
and will be the subject of a future publication. 

As a final remark, we point out that our formulation 
describes the particles as ideal chains. It is therefore 
natural to wonder how the results may be generalized to 
other species. A partial answer can be built following 
the "surface tension" argument presented below. It is 
known that the first order correction in Eq. H23|l for self- 
avoiding walks and spherical beads is comparable to that 
found for ideal chains [s^, so that it is tempting to ex- 
tend our conclusions concerning particle accumulation at 
the wall. This is not completely obvious however as the 
second-order coefficient /3 in the development Eq. H23|) 
does not follow this trend: (3 is positive for hard spheres. 
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negative for self-avoiding walks, and zero for ideal poly- 
mers. Nevertheless, the physical picture presented in this 
report is expected, at least qualitatively, to remain valid, 
and will await for further numerical and experimental 
verification. 
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In addition, verifies the boundary condition 

ZAr(r|| , u(r||)) = that we also expand as a series of 



=ZAr(r||,u(r||)) 
=Zjv(r||,0)+M(r||) 



dZ 



N 



^ 1 , X2 d^2N 



(A.4) 
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The problem is then solved recursively. The idea underly- 
ing this derivation is to recast the problem of modulated 
boundary conditions into a series of simpler problems in 
planar geometry. The calculations are performed up to 
second order, and we identify the corresponding condi- 
tions 



APPENDIX: PERTURBATIVE CALCULATION 

OF Zn 

The fundamental object in polymer statistic is the 
Green function or propagator G'Af({r||, z}, {rj|, z'}) of the 
chain. It represents the statistical weight of the random 
walk that starts at {r|| , z} and ends at {r|| , z'} in N steps. 
The Green function solves the Edwards equation |41(| 



on 



Gn = S{n)5{^ - r\)Siz ~ z') , (A.l) 



with b = a/\/6, a being the monomer size. The sur- 
face impenetrability is expressed as a boundary condi- 
tion Gjv({r||,'w(r||)},{r'|,z'}) — 0. For the flat, reference 
case, the Green function can be factorized 



G(^)({r||,z},{r:|,z'}) = (4^i?^ 



-3/2 



X < exp 



{z-z'y 



exp 



exp 
{z + z'Y 



mi 



(A.2) 



with Rg = V Nh^ the radius of gyration. 

The partition function of a chain with one ex- 
tremity fixed at {r||,z} is defined as 2^jv(r||,z) = 
/ (ir'|dz'G'Ar({r|| , z}, {rJi, z'}), the integral running over 
all the space available for the free extremity. For in- 
stance, the partition function of a chain with one cx- 



from a flat wall is Z 



(0) 



N 



duexp{—u^ ) is the error 



tremity fixed at distance 
eii{z/2Rg), where erf(2;) = 
function. 

The starting point of the derivation is to write the 

partition function as a series Zn — J2i ^ where Z^'' 
is of order [H El i^. The linearity of Eq. iTOT) 
then implies that each term satisfies an Edwards equation 
individually 



d_ 

dn 



(A.3) 



4'^ (r 11,0) = 



4'^ (r 11,0) 



-" I'll 

oz 

dz^^^ 



dz 



-"TV 



dz"^ 



(A.5) 



(A.6) 



The solutions of Eq. (jA.3|) can then be formally ex- 
pressed in terms of the Green function G*-") of the prob- 
lem il 



zW(r||,z) =6^^ dn j d^t 



(0) 



dz 



4L(rli,0). 



(A.7) 

This convolution product is conveniently written in 
Fourier-Laplace representation as 

9(^(0) 

i«(q,z) = 62__(q^^^0)i«(q,0) , (A.8) 



where the Laplace transform F of a function F/y is de- 
fined as F = /p dne~""Fn- It is then fairly simple to 
obtain the partition function up to second order 



1 



Z(i)(q, z) = -Uq-^exp --^/7+{q^))^ 
1 



(A.9) 



2(2) (q, z) = -^cxp [--yJTW 



d\' 



(A.IO) 



The set of Equations (|A.9|I . (|A.10p provides the exact 
expansion of the partition function of a Gaussian chain 
in the vicinity of a non-planar surface. The average value 
Eq. H1U|I of the density profile defined by Eq. H12|l follows 
immediately. 
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